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We study the finite-temperature thermodynamics of a unitary Fermi gas. The chemical poten- 
tial, energy density and entropy are given analytically with the quasi-linear approximation. The 
ground state energy agrees with previous theoretical and experimental results. Recently, the gen- 
eralized exclusion statistics is applied to the discussion of the finite-temperature unitary Fermi gas 
thermodynamics. A concrete comparison between the two different approaches is performed. Em- 
phasis is made on the behavior of the entropy per particle. In physics, the slope of entropy gives 
the information for the effective fermion mass m* /m in the low temperature strong degenerate re- 
gion. Compared with m* /m « 0.70 < 1 given in terms of the generalized exclusion statistics, our 
quasi- linear approximation determines m*/m ~ 1.11 > 1. 
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I. INTRODUCTION 

In recent years, the strongly interacting fermion 
physics becomes the focus of theoretical and experimen- 
tal attentionfl|. This is much attributed to the rapid 
progress of the atomic Fermi gas experiments. 

By tuning the external magnetic field, one can control 
the S-wave scattering length a or interaction strength 
between two atomic fermions. The crossover from 
Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein con- 
densation (BEC) can be realized by the so-called Fes- 
hbach resonance 0- At the resonance point, the scat- 
tering length can be singular with the existence of a 
zero-energy bound sate. Although the scattering length 
is singular, the scattering cross-section is saturated as 
a 47r/fc^ (with k being the relative momentum be- 
tween two atomic fermions) due to the unitary property 
limit. The divergent scattering fermion thermodynamics 
is referred to as the unitary Fermi gas thermodynamics 
in the literature Dealing with the strongly interacting 
matter is related with a variety of realistic many-body 
topics. 

Usually, the thermodynamics of dilute fermion system 
is determined by the two-body scattering length a, parti- 
cle number density n and temperature T. In the unitary 
limit with a = ±oo, the dynamical scattering limit should 
drop out in the thermodynamic quantities. At unitarity, 
the dynamical detail should not affect the thermodynam- 
ics; i.e., the unitary fermion system can manifest the uni- 
versal properties Q. 

Due to lack of any small expansion parameter, the uni- 
tary Fermi gas provides an intractable problem in sta- 
tistical physics. The fundamental issue is on the zero- 
temperature ground state energy. Based on the dimen- 
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sional analysis, the ground state energy should be propor- 
tional to that of the ideal Fermi gas with a universal con- 
stant ^ = 1 -|- /3, which excites many theoretical and ex- 
perimental efforts. The world average value of ^ is 0.42 — 
0.460, 0, 0, 0, Recently, we have attempted a quasi- 
linear approximation method to explore the strongly 
interacting limit fermion thermodynamics 9]. The ob- 
tained ground state energy or the universal constant 
^ = I is reasonably consistent with some theoretical or 
experimental investigations. 

Generally, the finite-temperature thermodynamics is 
as intriguing as the zero-temperature ground state en- 
ergy. There have been several Monte Carlo finite- 
temperature calculations of a unitary Fermi gas 10, 11]. 
In the strongly correlation unitary fermions, the nonlin- 
ear quantum fluctuations/correlations compete with dy- 
namical high order effects. In the weak degenerate Boltz- 
mann regime, the nonlinear correlations make the second 
order virial coefficient 02 vanish. To a great extent, the 
vanishing leading order quantum correction reflects the 
intermediate crossover characteristics of a unitary Fermi 
gas [a]. 

Can the intermediate characteristics be described in 
another way? In [l^, [l3| , the generalized exclusion statis- 
tics was developed to describe the anyon behavior in 
the low-dimensional strongly correlation quantum sys- 
tem. Physically, the behavior of a unitary Fermi gas is 
between Bose gas and Fermi gas[10|. Similarly, the be- 
havior of anyons is also between bosons and fermions. 
Can one use the anyons statistics to describe the inter- 
mediate unitary Fermi gas? Recently, the generalized 
exclusion statistics has been generalized to describe the 
unitary Fermi gas thermodynamics^^, ITsll. As a hypoth- 
esis, the priority is that the thermodynamics at finite- 
temperature can be investigated quantitatively. 

From the general viewpoint of statistical mechanics, 
calculating entropy is not a simple task. Either in clas- 
sical or quantum theory, the entropy describes how the 
microscopic states are counted properly. From the quan- 
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turn degenerate viewpoint, the low-temperature behav- 
ior of the entropy is a characteristic quantity. For ex- 
ample, according to the Landau theory for the strong 
correlation Fermi-Liquid, the slope of entropy per par- 
ticle versus temperature is related with the effective 
fcrmion mass m* /m. In physics, the dynamical param- 
eter m* jm is very important for the phase separation 
discussion for the as ymm etric fermion system with un- 
equal populations P^TTtI fist- Like the universal constant 
^ = 1 -|- /3, the effective fermion mass m* /m is an another 
universal constant for the BCS-BEC crossover thermo- 
dynamics. Obviously, the physics beyond the mean- field 
theory should be reasonably well understood. 

Unlike the ground state energy or the universal con- 
stant ^ with the world average value ^ « 0.44, the ef- 
fective fermion mass is an unknown parameter up to 
now. For example, the effective fermion mass is esti- 
mated to be m* jm k, 1.04 with a quantum Monte Carlo 
calculation 1 16]. A quantitative study of the phase dia- 
gram at zero temperature along the BCS-BEC crossover 
using fixed-node diffusion Monte Carlo simulations shows 
m* /m K, 1.09fl3|- A many-body variational wave func- 
tion with a T-matrix approximation leads to a larger 
value m*/m « LlTtlSj- What is the exact value of 
m* /ml 

In a quantitative way, we make a comparative study 
for the finite-temperature thermodynamic properties of 
the unitary fermion gas with the two formulations. The 
behavior of entropy per particle based on the quasi-linear 
approximation and the generalized exclusion statistics is 
discussed in detail. Indirectly, the effective fermion mass 
is determined from the entropy. The results are further 
compared with the Monte Carlo calculations. 

The paper is organized in the following way. In SecHIl 
the relevant thermodynamic expressions are given by the 
quasi-linear approximation. Correspondingly, the ther- 
modynamics given by the generalized exclusion statis- 
tics is presented in Sec lIIII The numerical calculations 
and concrete comparisons between the two methods are 
given in Sec lIVl In this section, the entropy per parti- 
cle and corresponding effective fermion mass m* /m are 
discussed. In Sec|Vl we present the conclusion remarks. 



tained grand thermodynamic potential ri(T, /x) or pres- 
sure P — —Q,/V can be described by the two coupled- 
parametric equations through the intermediate variable- 
effective chemical potential //* 



^ = Tj/5/2 ^ + ^^n' 

A'^ ' TO 
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In the above equations, A = \/2-k j (mT) is the ther- 
mal de Broglie wavelength and to is the bare fermion 
mass(with natural units ks — fi ^ i throughout the pa- 
per). 

The effective chemical potential /x* is introduced by 
the single-particle self-consistent equation, /i* makes 
the thermodynamic expressions appear as the standard 
Fermi integral formalism 



„v-l 



dx 



(3) 



where T(v) is the gamma function, and z = 1'^ is 
the effective fugacity. For example, the quasi-particle 
Fermi-Dirac distribution function gives the particle num- 
ber density according to 



(4) 



In the coupled equations Eq.(IT]) and Eq.®, the short- 
hand notations are defined as 



dn 



TT-: ■ (5) 



The shift term oc fir characterizes the high order non- 
linear contributions, which strictly ensures the energy- 
momentum conservation law. In the nonlinear approxi- 
mation, this significant high order correction term can be 
fixed in a thermodynamic way. It is worthy noting that 
the term oc /x^ can be exactly canceled by each other in 
the Helmholtz free energy density 



II. THERMODYNAMICAL QUANTITIES 
GIVEN BY STATISTICAL DYNAMICS WITH 
QUASI-LINEAR APPROXIMATION 

Strongly correlated matter under extreme conditions 
often requires the use of effective field theories in the 
description for the thermodynamic properties, indepen- 
dently of the energy scale under consideration. In the 
strongly interacting system, the central task is how to 
deal with the non-perturbative fluctuation and correla- 
tion effects. In Ref.|^, a quasi- linear approximation is 
taken to account for the non-local correlation effects on 
the unitary Fermi gas thermodynamics. 

With the quasi-linear approximation method, the ob- 



F 
V 
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where V is the system volume. However, the high order 
correlation term oc /x^ can be obtained in terms of the 
thermodynamic relations Q 



N 



and 



9£ 

9^, 



1^ ^ in 



dl 
dn 
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Comparing those obtained from Eq.([7]) and Eq.® with 
Eq.(IT]) and Eq.Q, the explicit expression of ji^ is 



dn 



2'Ka 



(9) 



The integrated expressions of the pressure and chemi- 
cal potential for the unitary Fermi gas are 

f3/2{z') , T/|/2(2')/-1/2(^') 



H = fi 



T 



fiMz') 



In the quasi-linear approximation, the auxiliary im- 
plicit variable fi* is introduced to characterize the non- 
linear fluctuation/correlation effects. As indicated by 
Eq. (fTO| and Eq.fTTjl. the ^* or z' makes the realistic 
grand thermodynamic potential Q{T,fi) appear as the 
set of highly non-linear parametric equations, which can 
be represented by the standard Fermi integral. By elim- 
inating the auxiliary variable the equation of state 
will uniquely be determined. 

From the underlying grand thermodynamic potential- 
partition function, one can derive the analytical expres- 
sions for the entropy density s = S/V and internal en- 
ergy density e = E/V. The following partial derivative 
formulae will be used 



dT 



dT 



/ J, 



4, 



\ dT J V ^^^* 
The entropy is derived according to 



s _ 1 /dP 
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Correspondingly, the explicit energy density expression 
is calibrated to be 



3T( , /|/2(^') /|/2(^')/-l/2(^T 

A3 ' 2A/2(z')+ 2/3/2(z') 



(14) 



Essentially, the entropy density includes the high or- 
der nonlinear contribution. What we want to emphasize 
is that the third law of thermodynamics is exactly en- 
sured as expected. The analytical analysis indicates that 
the energy density at zero-temperature gives the dimen- 
sionless universal coefficient according to ^ = fi/Ep = | 



or E/{^NEp) = ^, where the Fermi energy is Ep = 
(Sir^n)^^^ / (2m) and Tp is the Fermi characteristic tem- 
perature in the unit Boltzmann constant. The universal 
coefficient ^ = | has attracted much attention in the 
literature and is reasonably consistent with some Monte 
Carlo calculations [3. [lol| . 



III. THERMODYNAMICS GIVEN BY THE 
GENERALIZED EXCLUSION STATISTICS 

A. Generalized exclusion statistics 



T he g eneralized exclusion statistics is proposed in |12l | 
and [131 . If the dimensional of the Hilbert space is d and 
the particle number is N, then d and N are connected 
by Ad = —g A N, where the shift of the single-particle 
states number is Ad. The shift of the particle number 
for identical particle system is AA^ and g is a statistical 
parameter, which denotes the ability of one particle to 
exclude other particles in occupying single-particle state. 
When g = the intermediate statistics returns to the 
Bose-Einstein statistics and 5 = 1 to the Fermi-Dirac 
statistics. 

For anyons, the number of quantum states W oi N 
identical particles occupying a group of G states are de- 
termined by the interpolated statistical weights of the 
Bose-Einstein and Fermi-Dirac statistics. A simple for- 
mula with the generalized exclusion statistics is used to 
describe the microscopic quantum states [l3| 



W 



[G + (N - - g)]\ 
Nl[G - gN - {1 ~ g)]l- 



(15) 



One can divide the one-particle states into a large num- 
ber of cells with G ^ 1 states in each cell, and calculate 
the number with A^^ particles in the i-ih cell. The total 
energy and the total number of particles are fixed and 
given as 



(16) 



with Si defined as the energy of particle of species i. By 
generalizing Ea. (jl5p . we have 



w = l[ 



[G, + (Af,-l)(l-.g)]! 
A^,;![G.-5A^,-(l-5)]!' 



(17) 



We consider a grand canonical ensemble at tempera- 
ture T. For very large » 1 and A^^ ^ 1, using the 
Stirling formula InA^! = N{lnN — 1), and introducing 
the average occupation number defined by Ni = Ni/Gi, 
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one has 



In = ^ In 



[G, + (iV,-l)(l-g)]! 



N,\[G,-gN,-{\-g)]\_ 
- [^'' (1 + (1 - 9)N,) InG, (1 + (1 - .9)^0 

- G,(l - 5iV,) lnG,(l - gN.i) - G^N, InG^N,] . 

(18) 

Through the Lagrange multipUer method, the most 
probable distribution of Ni is determined by 



d 



(19) 



with chemical potential /i. It follows that 

N^e^^'-t')/T = [l + (l-g)7V,](i^9)(l-giV,)9. (20) 

Setting LOi — 1/ Ni — g, we have the anyon statistical 
distribution 



uj^+ g 

where uj obeys the relation 

One can define of w at e = with Eq. ([22] 
/. = -TlnK(l+LOo)^-^]. 



(21) 



(22) 



(23) 



The relation between /i and T has been established indi- 
rectly through ujQ and g. From Eq. (|22p . the uj and loq are 
related with each other through single-particle energy e 



e Tin 



LOQ 



1+UjQ 



which gives 



(24) 



(25) 



For T = 0, the average occupation number can be 
explicitly indicated as 

which is quite similar to the Fcrmi-Dirac statistics. 



B. Particle number and energy densities 

In the anyon statistics, the density of states is also 
given by 



D{e) = a(2m)3/Vei/V(47r2), 



(27) 



where a is the degree of the spin degeneracy and m is 
the bare fermion mass. 

At T = 0, the particle number is explicitly given by 



9 Jo 



(28) 



where Ep is related with the Fermi energy Ep through 
Ep = g'^/^Ep. With the Ep symbol, the system energy 
can be represented as 



E = - 
9 JO 



eD{e)de ^ -g^/^NEp. (29) 
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As we will see, once g is fixed, one can discuss 
the general finite-temperature thermodynamic proper- 
ties. Therefore, the essential task in the generalized ex- 
clusion statistics is fixing the statistical factor g. This 
can be determined by the zero-temperature ground state 
energy or the universal constant ^ according to ^ = g^^^. 
Various theoretical or experimental attempts have been 
made in the literature for determining the ground state 
energy. With the universal coefficient ^ — |[9], the ex- 
pected statistical factor can be identified to be g = ^. 

For the general finite-temperature scenario, the parti- 
cle number and energy can be rewritten as 



TV = 



E = 



' D{e)de 
uj + g ' 
eD{e)de 
uj + g 



(30) 
(31) 



By replacing Eq.l^-Eq.l^ and Eq.^ into Eq.iP]) 
and Eq. ([3T|) . one can have 



T 

2 V7> 
E 3 



3/2 



a(wo) = 1, 

\ 5/2 



NEp 
a(cjo) 



duj 



6(cjo) = 



duj 



In 



1 + W 
1 + CJo 



(32) 
(33) 

1/2 



1 + Wo 



-| 3/2 



Eq. (j32p determines ujq for a given temperature T. 
E/{NEp) can be obtained by a given uiq through Eq. ([55)) . 

For giving the explicit entropy density expression with 
the generalized exclusion statistics in the next subsection, 
let us make further discussion for the energy density. By 
eliminating N with Eq. ((28| and Eq. pS)) . the energy can 
be alternatively expressed as 



(34) 
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The partial derivative of the internal energy E to T for 
fixed fi is given by 



dT 



47r2 



(35) 



Furthermore, the variable ujo of the integral function 
&(a;o) can be converted into and T through Eq. 



b{ujo,fi,T) 



dio 



w(l + Lo) 



3/2 

(36) 



Therefore, one can have 



=^lnK(l + c.o)^-^]a(^o). (37) 



FIG. 1: The internal energy per particle versus the rescaled 
temperature. The solid curve denotes that for the ideal Fermi 
gas, and the short-dashed one is that given by the quasi-linear 
approximation. The long-dashed curve represents the result 
in terms of the generalized exclusion statistics model. The 
dots and solid squares are the Monte Carlo calculations 
and respectively. 



C. Entropy per particle 

Due to the scaling properties, the thermodynamics of 
a unitary Fermi gas also satisfies the ideal gas virial the- 
orem 



p = 



2E 
3V' 



(38) 



According to the thermodynamic relation for the en- 
tropy S and pressure P, one can have 



3 Ur 



(39) 



By substituting Eq.dSS]) and Eq.^ into Eq.(I3SI), the 
explicit expression for the entropy per particle is derived 
to be 



based on the quasi-linear approximation and the gener- 
alized exclusion statistics have similar analytical proper- 
ties; i.e., the internal energy increases with the increase 
of temperature. The two approaches both show that the 
energy density of a unitary Fermi gas is lower than that 
of the ideal Fermi gas. However, the shift of the internal 
energy given by the quasi-linear approximation is more 
quicker than that determined by the generalized exclu- 
sion statistics model. 




I = ^ (^)'^^(^o) +lnK(l +^o)^-^], (40) 
where luq is given by Eg. ([5^ for a given T. 



FIG. 2: Physical chemical potential versus the recaled tem- 
perature. The line-styles are similar to Figure. 1. 



IV. NUMERICAL RESULTS AND 
COMPARISONS 

Based on the above analytical expressions, we will give 
the numerical results. 



A. Internal energy and chemical potential 

From Eq. and Eq. ([55)1 , the energy per particle ver- 
sus the rescaled temperature can be solved. As indicated 
by FiglU the internal energies for the unitary Fermi gas 



With Eq.(l32]) and Eq.(l23l), we have also shown the 
chemical potential versus the rescaled temperature in 
Figl21 The chemical potential given by the two for- 
malisms decreases with the increase of temperature. The 
departure of them is getting bigger with the increasing 
temperature. 

The results for the energy per particle shown in Fig[T] 
and FigO in terms of the two different analytical ap- 
proaches are reasonably consistent with the Monte Carlo 
calculations El, while the chemical potential differs 
explicitly from the Monte Carlo result [ll| for T/Tp > 
0.8. 
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B. Entropy 
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FIG. 3: Entropy per particle versus the recaled temperature. 
The hne-styles are similar to Figure. 1. The Monte Carlo sim- 
ulation result is extracted from Ref. 10]. 

With Eg. ([5^ and Eq.fgD]), the entropy per particle 
curve versus the rescaled temperature is presented in 
FiglH The quasi-linear approximation predicts that the 
curve is higher than that of the ideal Fermi gas, while the 
generalized exclusion statistics model gives lower values 
compared with that of the ideal Fermi gas. With the 
increase of temperature, the entropy per particle given 
by the generalized exclusion statistics is getting closer to 
and almost overlaps with that of the ideal Fermi gas. In 
terms of the quasi-linear approximation, the ratio of en- 
tropy to that of the ideal Fermi gas approaches a constant 
in the Boltzmann regime. 

Especially, in the low-temperature strong degenerate 
regime, the slope of the entropy per particle versus the 
scaled temperature given by these two approaches is dif- 
ferent. The low-temperature behavior is determined by 
the effective fermion mass according to the Landau the- 
ory of strongly correlated Fermi-liquid. In turn, from 
the entropy curve, one can derive the effective fermion 
mass indirectly. The careful study shows that the quasi- 
linear approximation indicates m* /m « 1.11 > 1, while 
the latter predicts m* /m w 0.70 < 1. Compared with 



the latter, the quasi-linear approximation result is more 
consistent with the Monte Carlo calculations m* /m ~ 
1.04- 1.09 [H [13. 

V. CONCLUSION 

In terms of the quasi-linear approximation method and 
generalized exclusion statistics model, the internal ener- 
gies, chemical potentials and entropies of a unitary Fermi 
gas have been analyzed in detail. The two different ap- 
proximations give similar behavior for the internal ener- 
gies and chemical potentials of a unitary Fermi gas. 

The entropy is an important characteristic quantity 
in statistical mechanics. The entropy by the quasi- 
linear approximation is higher than that of the ideal non- 
interacting fermion gas. In the Boltzmann regime, the 
entropy curve given by the generalized exclusion statis- 
tics gets closer towards and almost overlaps with that 
of the ideal Fermi gas. The entropy given by the quasi- 
linear approximation is getting far away from that of the 
ideal Fermi gas and the ratio of entropy to that of the 
ideal Fermi gas approaches a constant. 

According to the quasi-particle viewpoint of the Lan- 
dau Fermi-Liquid theory, the slope of entropy per par- 
ticle determines the effective fermion mass in the low- 
temperature strong degenerate region. The numeri- 
cal analysis demonstrates that the generalized exclusion 
statistics model gives m* /m ~ 0.70 < 1. The developed 
quasi-linear approximation predicts m* /m ~ 1.11 > 1, 
which is closer to the updating Monte Carlo investiga- 
tions. 
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